Generation of macroscopic superposition states in ring superlattices 
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Ultracold bosons in rotating ring lattices have previously been shown to form macroscopic superpositions 
of different quasi-momentum states. We demonstrate that the generation of such kind of states using slightly 
non-uniform ring lattices has several advantages: the energy gap decreases less severely with the number of 
particles, the sensitivity to detunings from the critical rotation frequency is reduced, and the scheme is not 
limited to commensurate filling. We show that different quasi-momentum states can be distinguished in time- 
of-flight absorption imaging and propose to probe correlations via the many-body oscillations induced by a 
sudden change in the rotation frequency. 
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I. INTRODUCTION 

Macroscopic superposition states, to which we refer in this 
paper as (Schrodinger) cat states, are one class of strongly- 
correlated states which has attracted much interest since the 
early days of quantum theory [1|. Apart from foundational 
questions they have been shown to enable precision mea- 
surements at the Heisenberg limit [3| and serve as resources 
for various quantum information tasks (4). Efficient cat state 
production is therefore an issue of theoretical and practical 
importance. 

Among theoretical proposals for creating cat states with 
Bose-Einstein condensates |5j [6] [71 [8) one is to load ultra- 
cold bosons into rotating ring lattices |9] \W[ QT) . At a critical 
rotation frequency the two lowest-lying states are the sym- 
metric and anti-symmetric superpositions of different flow 
around the loop I91 ITT1 . To observe these strongly-correlated 
states is however largely impossible as their energy differ- 
ence decreases exponentially in the number of particles ifTTl . 
On the other hand, superpositions of macroscopically distinct 
flow have been observed in experiments with superconducting 
quantum interference devices (SQUID) lfT2l [131 fl4l . i.e. su- 
perconducting rings with the superconductivity deliberately 
destroyed at one or several places in the loop fl5l . 

In this paper we address the question whether cat state pro- 
duction with ultracold bosons in rotating ring lattices can be 
improved by varying the coupling due to tunneling along the 
ring. After introducing the twisted Bose-Hubbard Hamilto- 
nian for this problem, we diagonalize analytically the single- 
particle Hamiltonian for ring superlattices where the tunnel- 
ing matrix elements alternate along the ring. We show that 
in slightly non-uniform rings with weak on-site interactions 
there is a critical rotation frequency at which, similar to the 
case of a uniform ring lattice, the ground and first excited 
states become cat-like superpositions of the two lowest-lying 
degenerate single-particle modes. We demonstrate that their 
energy difference scales less severely with the number of par- 
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tides, the superposition states are less sensitive to detunings 
away from the critical rotation frequency compared to the uni- 
form case, and cat state production is not limited to com- 
mensurate filling of the lattice. Finally, we discuss how one 
can probe cat-like correlations via many-body oscillations in- 
duced by a sudden change in the rotation frequency and show 
that different quasi-momentum states can be distinguished in 
time-of-flight absorption images. 



II. THE HAMILTONIAN 

We consider a system of N ultracold bosons with mass M 
confined in a ID ring superlattice of an even number of sites 
L with lattice constant d. The ring is rotated about its axis (z 
axis) with angular velocity il. A theoretical proposal how to 
experimentally achieve such a ring-shaped optical lattice can 
be found in Ref. Ifl6l . In the rotating frame of the ring the 
many-body Hamiltonian is given by ifTTlfTTll 
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where a is the s-wave scattering length, V^x) the lattice po- 
tential, L z the angular momentum, and x the 3D spatial co- 
ordinate vector. As the lattice potential V'(x) is a ID ring 
superlattice in the x — y plane, it confines the motion along 
the z axis as well as the radial direction in the x — y plane so 
strongly that only the motion of the atoms along the ring has 
to be taken into account. 

Assuming that the lattice is deep enough to restrict tun- 
neling between nearest-neighbor sites and the band gap 
is larger than the rotational energy, the bosonic field op- 
erator <3E> can be expanded in Wannier orbitals confined 
to the first band $ = J2j^] W j( x ) with W j( x ) = 
exp \=f±- /* A(x') • dx'] Wj(x). Here Wj(x) are the Wan- 

nier orbitals of the stationary lattice centered at the site j, 
A(x) = Q,z x x an effective vector potential and dj the 
bosonic annihilation operator of a particle at site j. In terms of 
these quantities, the many-body Hamiltonian can be written, 
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j=o I otherwise 

( 4 ) 

we obtain the single-particle Hamiltonian H s? in the quasi- 
momentum basis, i.e. 
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FIG. 1 : (Color online) One-dimensional ring lattice with L = 4 sites 
and alternating tunneling matrix elements J and t. In the rotating 
frame an effective phase twist 9 is induced. 
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Here the two-by-two matrices are given by 

/ -(J + t)cos(e-^) -i(J-t) S in(9-^l) 
q \ +i(J — t) sin [9 — +( J + t) cos (0 - ^a) 

(6) 

and represent the coupling between the quasi-momentum 
states. We can diagonalize the single -particle Hamiltonian |5]) 
via a unitary basis transformation 
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up to on-site diagonal terms which we neglect for simplicity, 
as 10 DH 03 
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Here 
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is the number operator at site j, 9 
is the effective phase twist induced by the gauge field, 
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Jj is the hopping en- 



ergy between nearest neighbor sites j and j + 1: J, 



2M 



^(X) 



and J7 the on-site inter- 



action energy: U = 



— 4Trah 



M 



The case of a uniform ring, i.e. Jj = J for all j, has been 
discussed before. In this paper we generalize to non-uniform 
ring lattices where we will focus on superlattice structures 
where Jj = J for j even and Jj — t for j odd (see Fig. [TJ. 
In this case the single-particle problem can be solved analyti- 
cally. 



III. THE NON-INTERACTING SYSTEM 

We shall first study the properties of the Hamiltonian Q in 
the limit of vanishing on-site interaction U = 0. 
Using the decompositions 
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where the single-particle energies are given by 



Ef = ±\\J 2 + t 2 + 2 Jt cos ( 29 - 



which for t = J simplify to ifTTI 



£7* = ±2J 
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(10) 



(11) 



In the uniform ring, t/J = 1, the eigenstates of the single- 
particle Hamiltonian H sp are quasi-momentum states. At cer- 
tain phase twists 9 they are twice degenerate. For example, 
for L = 4 sites the quasi-momentum states \q = 1) and 
|<7 = —1) are degenerate at 9 = 0, whereas at 9 = ir/4 the 
states \q = 0) and |g = 1) as well as \q = 2) and |g = —1) are 
degenerate. Reducing the symmetry of the ring by choosing 
t J the quasi-momentum states which differ by L/2 quasi- 
momentum units are coupled by the single-particle Hamil- 
tonian d5|, so that the quasi-momentum states \q = 1) and 
\q = —1) hybridize and the degeneracy at 9 = is lifted, 
i.e. E+ - El = 2(J- t). 

At 9 — ir/4 however the degenerate quasi-momentum 
states are not coupled by the single-particle Hamiltonian ([5J, 
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FIG. 2: (Color online) Single-particle spectrum for L = 4 sites with 
t/J — 1 (left) and t/J = 0.7 (right) as a function of the phase twist 
8. In both cases there are level crossings at 9 = 7r/4. 



so that the degeneracy is present also in the non-uniform case. 
This remains true for arbitrary L, i.e. Eq = -E£/ 4 at 8 = tt/4, 
but for L 4 these states are not the ground states of the sys- 
tem. In Fig. [2] we plot the single-particle spectrum for L = 4 
sites as a function of the effective phase twist 8. It shows level 
crossings at 8 = tt/4 both for t = J as well as t ^ J. In the 
following we will refer to 8 = tt/4 as the critical phase twist, 
since - as we will demonstrate below - weak on-site interac- 
tions lift the degeneracy at 9 = tt/4 and lead to the formation 
of strongly-correlated states in the many-body system. 



IV. EFFECTIVE MANY-BODY HAMILTONIAN 



In this section we derive an effective many-body Hamilto- 
nian to study the properties of the many-body system in the 
limit of weak on-site interactions. In the absence of interac- 
tions U = the ground state of a bosonic many-body sys- 
tem is the state with all bosons occupying the lowest-energy 
single-particle state. At the critical phase twist 8 = tt/4 the 
single-particle Hamiltonian <|3j has, however, two degenerate 
single-particle states, so that there is a N + 1 -dimensional 
degenerate subspace at iV-particle level. A convenient basis 
for this subspace is given by the Fock states |n, N — n) with 
< n < N, where n particles are in the single-particle state 
of energy Eq and N — n particles in the one of energy E^u, 

respectively. For weak interactions NU/L <C 2\/ J 2 + t 2 this 
subspace is the low-energy sector of the many-body problem 
for all phase twists 8 and tunneling strength ratios t/J. 



Starting from the interaction Hamiltonian 
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we first restrict the sum in Eq. ( |T2] > to the relevant modes 
ki e {0,L/4,L/2,3L/4} with i = {1,2,3} and second use 
the inverse of Eq. ^ to express the operators b q in terms of 
the operators c q . Keeping only terms within the low -energy 
subspace we obtain the effective Hamiltonian up to first order 
in the on-site interaction strength U 



H, 



eff 



(E h + E L/i h L/ ^j + — (2n Q h L/4 + N 2 -N) 
+ (^j^clclc L/i c L/i + H.c\ (13) 

where h q = c q c q are the number operators and the parameter 

rj = 2 (cos Qfo sm a o cos 2«l /a 

+ cos a L / 4 (— sin a£/ 4 )cos 2ao) (14) 

satisfies < r/ < 1 and simplifies for 8 = tt/4 to 

J 2 -t 2 
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The first bracket of Eq. ( 13 1 contains the contributions from 



the single-particle Hamiltonian ([9]), whereas the terms in the 
second and third brackets arise from the on-site interaction 
( 12 1. At the critical phase twist 8 = tt/4 the former are an 
unimportant zero-energy offset, whereas the terms in the sec- 
ond bracket shift the energies of the states in the subspace dif- 
ferently, e.g. they lead to an energy difference of U(N — 1)/L 
between the states |iV, 0) and \N — 1,1), while the states 
1 71, N — n) and \N — n, n) remain pairwise degenerate. The 
terms in the third bracket are off-diagonal in the Fock basis 
of the subspace and describe two-particle scattering between 
the two single-particle modes. As we will see in the next sec- 
tion they lift the remaining pairwise degeneracies in the many- 
body spectrum. In Fig. [3]we plot the many-body spectrum for 
L = N = 4, t/ J = 0.7 and U/J = 0.5 as a function of the 
phase twist 8 obtained from exact diagonalization of (|2]l and 
the effective Hamiltonian ( 13 I. We see that for these parame- 



ters the effective Hamiltonian ( 13 1 is an excellent approxima- 
tion for the low-energy sector of the many-body Hamiltonian 



A. Cat-like superpositions in the limit i/Jftil 

Let us now determine the ground and first excited state for 
slightly non-uniform rings !/J« 1, close to the critical phase 
twist 8 w 7r/4. Since the terms in the second bracket of 



Eq. ( 13 i increase the energy for all states in the subspace apart 
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FIG. 3: (Color online) Many-body spectrum for L — N — 4, t/J = 
0.7 and [//J = 0.5 as a function of the phase twist 8 obtained from 
exact diagonalization (dotted lines) of the many-body Hamiltonian 
Q as well as from the effective Hamiltonian ([13} (solid lines). 



from \N, 0) and |0, N) and the coupling between the states is 
weak (as the coupling r\ is small in this limit), we project the 
effective Hamiltonian ( 



13 i onto the subspace spanned by these 



two nearly-degenerate lowest-energy states ifTTIl . 

As there is no direct coupling between \N, 0) and |0, N) we 
calculate the total coupling through intermediate states using 



perturbation theory. After eliminating the intermediate states 
we obtain the following 2x2 Hamiltonian 
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where AE is the energy difference between the states \N, 0) 
and |0, N) caused by the detuning of the phase twist from 
resonance A8 = 9 — it /A, i.e. 



AE = N(E- /i -Ev) 
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and A is the coupling between the states \N, 0) and |0, N) due 
to the off-diagonal terms of the effective Hamiltonian (13) . 
As the latter only directly couples the states \n, N — n) and 
\n ± 2, N — n =p 2), the first non-vanishing order is given by 



A = 
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(18) 

Here, Ej = jfc (2j(N - j) + N 2 - N) is the diagonal in- 
teraction energy shift. Finally, the term £ in Eq. ( 16 1 is the 
energy shift induced by the non-diagonal terms in the effec- 
tive Hamiltonian ( p"3j ) 
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As £ is just a constant in the 2x2 subspace we do not need to 
evaluate it explicitly and neglect it hereafter. 



The ground state of the two-by-two Hamiltonian (16i is 
given by 



c \N,0)+i N / 2 c N \0,N) 
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where the ratio of its amplitudes is given by 
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and the energies of the ground and excited states are E± = 
±^AE 2 + |2A| 2 /2. We see that to obtain a cat-like super- 
position, i.e. co/cjv ~ 1, the energy difference must not dom- 
inate over the coupling \AE\ <C |2A|. In this limit the energy 
gap E + - E_ is given by E + - E_ w |2A|. 

In Fig. |4] we plot the energy gap between the two lowest- 
energy states of the many-body spectrum as a function of the 
on-site interaction U/J for L = N = 4, t/J = 0.7 and 
9 = 7r/4. We find that the analytic formula ( 18 i agrees very 



and at small interaction strengths U < 0.2 also with the exact 
diagonalization of the full many-body Hamiltonian 



Let us compare our result ( 18 1 with the coupling A in the 
case of a uniform ring lattice (t = J). Unfortunately, an an- 
alytic expression equivalent to Eq. (TE\ has not been found 



well with the diagonalization of the effective Hamiltonian ( 13 i 



in this case but it is known ifTTI that the coupling A in first 
non-vanishing order for unity filling N/L = 1 is propor- 
tional to A/J oc {jfjj) N ~ l - In both cases the coupling A 
decreases exponentially with the number of particles N. This 
is because multiple two-particle scattering processes are the 
microscopic origin of the coupling A and transitions between 
the two configurations \N, 0) and |0, N) are thus highly off- 
resonant N/2-th and N — 1-th order processes, respectively. 
Consequently, the energy gap vanishes in the thermodynamic 
limit N — > oo and cat state production both in uniform and 
slightly non-uniform ring lattices is restricted to modest num- 
bers of atoms. However, since the energy gap in first non- 
vanishing order is proportional to U (r]/2) N / 2 for slightly non- 
uniform rings as compared to U(U / J) w_1 for uniform rings, 
it can be one order of magnitude bigger in the relevant param- 
eter regime. Furthermore, we point out that the appearance of 
cat states in slightly non-uniform rings is not limited to com- 
mensurate filling as long as the number of particles N is even. 
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FIG. 4: (Color online) Energy gap for L = N = 4, t/J = 0.7, 
6 = 7r/4 as a function of U/ J from exact diagonalization (blue 
points), effective Hamiltonian (red line) and the gap formula (yellow 
line). 
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FIG. 5: (Color online) Energy gap as a function of particle number 
N with t/J = 0.7 (upper curve) and t/J — 1 (lower curve) for the 
ring of L — 4 sites and U/J — 0.5. The solid lines are exponential 
fits to guide the eye. 



In Fig. [3] we plot the energy gap L = 4 and U/J = 0.5 
as a function of number of particles N for t/J = 0.7 and 
t/J = 1. The solid line is an exponential fit to the points 
to guide the eye. We see that by changing from t/J = 1 
to t/J = 0.7 the energy gap increases by over one order of 
magnitude for numbers of particles N > 8. As experiments 
are typically limited to lsec and typical hopping energies are 
of order of 0.05i? r ~ 10 2 Hz, the detection of cat states by 
non-equilibrium dynamics is limited to A/J = 10~ 2 . This in 
turn restricts our cat-state production scheme to a few tens of 
particles in contrast to less than ten particles for uniform ring 
lattices (see Fig. [5]). 



B. Product states in the limit t/J — 

Let us now try to understand the properties of the slightly 
non-uniform ring by looking at the opposite limit t/J=0 



where the system becomes an array of L/2 isolated double- 
well potentials. In this section we will focus on the case of 
L = 4 sites for notational convenience, but our arguments are 
valid for all L and N. 

Since the double-wells are isolated for t/J = 0, we can 
treat the atoms in each double-well separately. Without inter- 
actions U = the states \n, N — n) (0 < n < N) with n 
atoms in the ground state of the first and N — n in the ground 
state of the second double-well have the same single-particle 
energy. In first-order perturbation theory the on-site interac- 
tion shifts the ground state energy of a double-well with n 
particles by U/2 ■ n(n — l)/2 and thus breaks this degeneracy. 
The energies of the states \n, N — n) become 
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(22) 



As expected, repulsive interactions favor states of equal num- 
bers in the two double-wells, so that the states | N/2, N/2) and 
\N/2 ± l,N/2 =p 1) are the lowest-lying states. The energy 
difference between the ground and the two-fold degenerate 
excited state is thus U/2. We note that the opposite sign of 
the energy shift due to the on-site interaction in Eqs. ( 13 1 and 



( 22 1 stem from the two different single-particle basis sets used 
to span the two-mode Fock states EUll . 

In Fig. [6] we compare the excitation energies of the effective 
Hamiltonian ( p"3j ) with the full many-body Hamiltonian |2]) as 
a function of t / J and find that they agree for all ratios t/J. 
Furthermore, we see that in the limit t/J — > the excitation 
energies agree with the picture of the isolated double-wells. 
In the inset of Fig. [6] we plot the overlap of the ground state 
of the system with the ground state of the isolated double- 
wells t/J = as well as the uniform ring t/J = 1. We find 
that as the ratio t/J decreases, the overlap with the cat-like 
superposition of two quasi-momentum states decreases while 
the overlap with the ground state of the isolated double-wells 
increases. 

We further note that in the limit t/J — > the physics is in- 
dependent of the number of particles N and the applied phase 
twist 9 (see Eq. ( 10 ). As the exact ground state of the ef- 
fective Hamiltonian ( 1 3 1 for arbitrary t/J and N is very non- 
trivial, a legitimate way to gain some physical insight about its 
properties is to interpolate between the two limits discussed 
above. Using this perspective one can attribute the less se- 
vere scaling of the energy gap with the number of particles, to 
hybridization of the cat and product states. 



V. BIGGER ENERGY GAPS AND WEAKER SENSITIVITY 

A fair comparison between cat state production in uniform 
and non-uniform ring lattices has to contrast the gain in energy 
gap with the degradation of the cat fidelity or "cattiness". To 
quantify the "cattiness" of a quantum state we consider two 
measures: (i) we define the cattiness of the first kind as the 
overlap of the ground state with the quasi-momentum cat state 
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FIG. 6: (Color online) Excitation energies for L — 4, N = 16 
and 8 = ir/4 as a function of t/J from the effective Hamiltonian 
l|T3} (solid line) and the full Hamiltonian (|2j (dotted line). The in- 
set shows for the same parameters the overlap of the ground state 
with the ground state of the isolated double- wells t/J = (red dot- 
ted line) as well as the uniform ring t/J=l (blue solid line) as a 
function of t/J. 



and (ii) we define the cattiness of the second kind as the over- 
lap of the ground state with the cat-like superposition of the 
two single-particle modes 
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which depends on the ratio t/J as defined in Eq. Q. 

In this section we demonstrate that while there is a trade- 
off between bigger energy gaps and cattiness of the first kind, 
non-uniform ring lattices offer bigger energy gaps at fixed cat- 
tiness of the second kind. Consequently, they are less sensitive 
to detunings of the rotation frequency which is an important 
practical advantage. On the other hand it is clear that to detect 
a cat state in two single-particle modes which do not have well 
defined quasi-momentum is more demanding experimentally. 

Both cat-state measures can be calculated perturbatively in 
the limit of weakly non-uniform rings t/ J w 1 and weak on- 
site interactions U/J <C 1. The quasi-momentum cat state 
|?/>i) is the ground state for t/J — 1 and U/J — > 0. Devi- 
ations are due to the off-diagonal part of the single-particle 
Hamiltonian in quasi-momentum representation Q as well as 
12 1. We obtain for the overlap 



the on-site interactions 



,; , , , v,2 N ( J -t 

K^>l =1 -T {iTt 



'Uy/N(N-1) 

2LVJ 2 + &~ 



(25) 



where |?/>o) denotes the ground state of the system at finite t/J 
and U I J. In case of the cattiness of the second kind, depletion 
is due to the off-diagonal terms of the effective Hamiltonian 



(13 1 as well as the part of the interaction Hamiltonian (12i 
which couples states within the subspace to states outside the 
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FIG. 7: (Color online) Energy gap versus cat-state measure of the 
first kind (dotted points) as well as of the second kind (solid lines) 
for t/J = 0.7 (upper blue) and t/J — 1 (lower red). The other 
parameters are L — 4, N = 16 and 9 — n/4. 



subspace of the effective Hamiltonian. We obtain 
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In Fig. [7] we plot the energy gap versus the two cat-state 
measures for t/ J = 0.7 and t/J = 1, respectively. The other 
parameters are L = 4, N = 16, 9 = 7r/4. We find that: (i) the 
larger U / J, the bigger the energy gap but the smaller are both 
cat-state measures. This is readily understood from Eqs. (18), 



( 25 1 and ( 26 1: stronger on-site interactions increase the cou- 
pling |A| and, therefore, increase the energy gap. However 
they also increase coupling to other states outside of the sub- 
space of i?2x2 an d the cat state is gradually depleted ifTOl . (ii) 
The smaller t/J, the bigger the energy gap but the smaller is 
the cattiness of the first kind at fixed U/J. This is in agree- 
ment with Eqs. ( 18 1 and (26 1. Physically, this means that the 



off-diagonal part of the single-particle Hamiltonian increases 
the coupling |A| but also dilutes the many-body correlation in 
the quasi-momentum basis. Consequentially, there is a trade- 
off between energy gap and cattiness of the first kind. We note 
that it might be a sizeable advantage to sacrifice part of the cat- 
like correlation (and therefore a weaker signal in the many- 
body oscillations we discuss in the next section) but make it 
visible on experimentally accessible time scales. Finally, we 
find that (iii) one can increase the energy gap by about an or- 
der of magnitude (for N = 16 particles) at fixed cattiness of 
the second kind when the ratio t/J is changed from t/J=l 
to t/ J — 0.7. While we believe that this is the natural measure 
of comparison between uniform and non-uniform ring lattices, 
we admit that probing the cat state in the single-particle basis 
is a non-trivial task, since - as we will show below - time- 
of-flight imaging maps more closely on the quasi-momentum 
basis. 

Increasing the coupling | A | with non-uniform ring lattices 
has important implications for the sensitivity to detunings of 
the rotation frequency away from the critical phase twist. In 
Fig. [8] we plot the cat-state measure of the second kind as a 
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FIG. 8: (Color online) Cat-state measure of the second kind as a 
function of t/J and the detuning AO while L = 4, N = 16 and 
U = 0.2 are fixed. 



FIG. 9: (Color online) Many-body dynamics |6o(i)| 2 and |6jv(t)| 2 
for L — 4 and N — 8 following a sudden change in phase twist 
from 9 = to 6 = tt/4 at (7/ J = 1 and t/J = 0.7. 



function of t/J and the detuning AO while L = 4, N = 16 
and U — 0.2 are fixed. We find that the system is much less 
sensitive to detunings away from the critical phase twist for 
t/J< 1, i.e. the requirements on the phase control precision 
AO are relaxed. To understand how this comes about let us go 



of the probability to be in the state \N, 0) and |0, N) is 



back to the two-by-two Hamiltonian ( 16 1. Similar to the two 



state Hamiltonian for a single particle in a double-well poten- 
tial, its ground state depends on the ratio of energy difference 
AE to the coupling energy |A|: for \AE\ 3> |A| the system 
is in one of the states \N, 0) (in the double-well analogy on the 
left-hand side) or |0, N) (on the right hand side) whereas it is 
in their superposition state in the limit where the coupling en- 
ergy dominates \AE\ -C |A|. As the coupling |A| in Eq. ( 18 i 
decreases exponentially with increasing number of particles 
TV, more control over the rotation frequency is needed to tune 
systems with more particles into resonance, i.e. to make the 
energy difference \AE\ between the states \N, 0) and |0,iV) 
smaller than their coupling |A|. This was identified as one of 
the barriers for creating large cat-like states ifTOll . As the gap 
increases when t/J is decreased (see Eq. ( 18 1), we conclude 



that using non-uniform ring lattices t ^ J can improve the 
situation on this issue. 



|6 (t)| 2 = |(iV,0|^))| 



1 



cos vt 
2 



(27) 



and \b N (t)\ 2 w 1 — |&o(i)| 2 where hv = E + — E_ denotes 
the energy gap and where the approximate signs indicate that 
we neglect the depletion due to the non-diagonal part of the 
single-particle Hamiltonian |5]l as well as the on-site interac- 
tion Hamiltonian ( p"2| ). 

In Fig. p\ we show the many-body dynamics, i.e. |6o(t)| 2 
and \bff(tj\ 2 for L = 4 and N = 8 following a sudden 
change in phase twist from # = 0to6> = 7r/4at time t = 
with U/J = 1 and t/J = 0.7. We see that the many-body 
oscillations following this sudden change occur at the fre- 
quency of the gap hv. These oscillations are modulated by 
an oscillation with frequency U(N — 1)/L which is the en- 
ergy difference between the nearly-degenerate ground states 
and the next lowest-lying intermediate states \2,N — 2) and 
\N — 2, 2). The amplitude of the oscillations is smaller than 
one, i.e. |&o(*)| 2 + l^ ) A r ( f )| 2 < 1> due to the depletion of the 
quasi-momentum cat state |(-0i|?/>o)| 2 ( see Eq. (25 i). 



VI. DYNAMICAL DETECTION OF CAT-LIKE 
CORRELATIONS 

In this section we propose to induce many-body oscillations 
by suddenly changing the applied phase twist 8 in order to de- 
tect the coherent superposition of two quasi-momentum states 
|JV,0> = (SjHvacVA/M and \0,N) = (S[ /4 ) w |vac)/vOV! 
at the anti-crossing of the many-body spectrum. 

If we assume the system is initially in the ground state of 
the full Hamiltonian at 9 = which will predominantly be 
the state |iV, 0), i.e. \tp(t = 0)) w \N,Q). Then, the phase 
twist is changed to 9 — ir/4, where the eigenstates are ap- 
proximately |±) w (\N, 0) ± |0, N))/V2. The time evolution 



VII. DETECTION 

In this section we study the spatial density profile of an 
atomic cloud when it is let to freely expand after turning off 
the ring lattice potential. Our analysis shows that states with 
different quasi-momentum states exhibit a distinctive time-of- 
flight pattern and that therefore the latter can be used as an 
experimental probe for detecting the many-body oscillations 
between the |0, N) and \N, 0) discussed in previous section. 

For this discussion we will consider the expansion of the 
wavefunction of a single atom initially confined in a ID lattice 
ring geometry. For simplicity we approximate the axial, radial 
and angular Wannier functions as Gaussians with width er 2 , ay 
and Pq(Jq, respectively, and localized in a ring of radius po at 



8 



ip(r,t = 0) = Ae 



(r-PoY 



L 



C n € 



(28) 



where A is a normalization constant. If at time t = the atom 
is in an eigenstate of the uniform ring with quasi-momentum 
h27rq /La then c„ = * e^ 2 ™*)/^. 

At sufficiently large times t ^> Mpo / h one can take the 
far-field limit and approximate the density of the expanded 
cloud as being proportional to the momentum distribution of 
the initial state, i.e. 



|V(r,^oo)| 2 



(29) 



where Q(r) = Mv/ht and |0(k)| 2 is the momentum distri- 
bution of the initial state. After integrating along the z axis, 
one can write a simple analytical expression of the transverse 
spatial density distribution in the limit of strongly confined 
Wannier orbitals (i.e. oy — > and a® — ► 0) 



\ip(p,Q,t -xx)\ = / \4>(r,t -> oo)\ 2 dz 



fMp 



-iQ(p) Po cos(2^a-e)) 



(30) 



V Lht 

In the limit of a large number of lattice sites L, the sum in 



Eq. ( 30 1 can be approximated by an integral and the time-of- 



flight profile reduces to 

l^( p ,e,;->oo)| 2 -> (^-J \J qa [Q(p) P o]\ 2 (3D 

where J qo (x) are Bessel functions of the first kind. This case 
corresponds to the rotationally symmetric case discussed in 
Ref. Ifl8l where the position of the zeros and maxima of the 
Bessel functions provide a full characterization of the density 
profile. A state initially with zero quasi-momentum for exam- 
ple will exhibit an interference peak at the origin while a state 
with non-zero quasi-momentum will exhibit a central hole. As 
the position of the first maximum of J qo [Q(p)po] is an increas- 
ing function of qo, the larger the initial quasi-momentum the 
wider the central hole. For a finite number of lattice sites the 
sum does not correspond exactly to a Bessel function and the 
momentum distribution does not become fully radially sym- 
metric. However, there is still a unique correspondence be- 
tween the position of the peaks in the absorption image and 



the initial quasi-momentum distribution (see also Ref. Ifl9l ). 
Consequently, the latter provides a means to experimentally 
determine the quasi-momentum of the wavefunction before 
the release. 



In Fig. 10 we show the numerically calculated time-of- 
flight images of different ring lattice geometries for two states 
with initial quasi-momentum qo = and L/4, respectively. 
These states correspond to the two types of initial distributions 
that one has to experimentally distinguish in order to probe the 
many-body oscillations we have discussed in the previous sec- 
tion. The figure shows the distinctive interference pattern of 
each of the two initial states and the asymptotic approach to 
radial symmetry as the number of lattice sites L is increased. 



VIII. CONCLUSION 

We have investigated cat state production with ultracold 
bosons in rotating ring superlattices and found the improve- 
ment compared to uniform ring lattices to be threefold: the 
energy gap between the cat-like ground and first excited states 
scales more favorably with the number of particles and the 
constraints on phase twist control and commensurate filling 
are relaxed. Finally, we have shown that the different quasi- 
momentum states can be distinguished in time-of-flight ab- 
sorption images and proposed to probe the cat-like correla- 
tions via the many-body oscillations induced by a sudden 
change in the rotation frequency. Since the exponential scal- 
ing of the energy gap with the number of particles remains, 
we suspect that cat state production in ultracold atomic gases 
is limited to modest number of atoms at least for systems 
with contact interactions. One possible way-out is to consider 
systems with long-range interactions such as polar molecules 
with dipolar interactions. 
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